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Abstract 

We employ the Sagdeev pseudo-potential method to investigate the propagation of nonlinear ion 
waves in a relativistically degenerate electron-ion plasmas. The matching criteria for existence of 
such nonlinear excitations are numerically investigated in terms of the relativity measure (relativis- 
tic degeneracy parameter) of electrons and the allowed Mach-number range for propagation of such 
waves is evaluated. It is shown that the electron relativistic degeneracy parameter has significant 
effects on nonlinear wave dynamics in superdense degenerate plasmas such as that encountered in 
white dwarfs and the cores of massive planets. 
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I. INTRODUCTION 



Recently, investigation of the plasma under extreme conditions, such as high temperature 
and pressure, has attracted outstanding attention because of numerous applications in labo- 



ratory 



as well as astrophysical environments js-O]. Under such conditions the quantum 
and relativistic effects are unavoidable and one must take into account the relativistic de- 
generacy pressure which arises due to the Pauli-blocking mechanism. The quantum effects 
become more important when the inter-fermion distances are much lower than the charac- 
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Chandrasekhar 1^, using 



teristic de Broglie thermal wavelength, h/ {2TimkBTy^'^ 
the Fermi-Dirac statistics for electrons, has developed a rigorous mathematical criteria un- 
der which a white dwarf can be considered as completely degenerate, i.e. zero-temperature 
Fermi-gas. The complete degeneracy may seem a simplifying assumption, however, it hold 
for many dense astrophysical entities such as white dwarfs jisl . [l^ . For a completely degen- 
erate plasma one can assume well defined Fermi quantities such as that for typical metals. 

The electron degeneracy in massive white dwarfs which opposes the catastrophic grav- 
itational inward pressure can give rise to relativistic degeneracy at some condition. Many 
studies in the framework of quantum hydrodynamics dealing with astrophysical models con- 
sider only the non-relativistic degeneracy cases 
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-2d]. However, in a typical white dwarf 
where the electron density can be as high as lO^^cm"^ the degeneracy can be relativistic 
and both quantum and relativity considerations should be taken into account. Semiclassical 
approach shows that the dynamics of nonlinear ion-acoustic propagations in an electron- 
positron-ion can be much different under relativistic degeneracy pressure from those under 



non-relativistic one 



23| . The presented investigation based on the quantum hydrodynam- 



ics model does not cover all range of the relativity parameter for plasma under consideration 
and are relevant only to either non-relativistic or ultra-relativistic cases. 

The aim of current investigation is to study the propagation of nonlinear excitations, 
namely, ion solitary waves (ISWs) and ion periodic waves (IPWs) with arbitrary amplitudes 
in relativistically degenerate zero-temperature electron-ion plasma. Although, it has been 
confirmed that the model of zero-temperature Fermi-gas is quite appropriate for most white 
dwarfs [12|, however, the full determination of large-amplitude propagation of ISWs and 



IPWs in such plasma with extreme temperatures requires the finite-temperature Fermi- 
Dirac model which accounts for the electron temperature that may be our future work. 
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The presentation of the article is as follows. The basic normalized plasma equations are 
introduced in section [III Nonlinear arbitrary-amplitude solutions are derived in section IIIII 
Numerical analysis and discussions is given in section IIVI and final remarks are presented in 
section |Vl 



II. BASIC HYDRODYNAMICS EQUATIONS 

We consider a completely degenerate dense plasma consisting of relativistic cold-electrons 
and dynamic cold-ions. It is noted that in degenerate plasmas the rate of electron-ion 
collisions are limited due to Pauli-blocking mechanism which allows only degenerate particles 
with energies limited to a narrow range around the Fermi-energies to interact, hence, the 
plasma may be considered almost collision-less. We use the conventional fluid hydrodynamics 
equations, which in dimensional form can be written as 

drij I driiVj _ n 
dt dx ^' 

dt dx ^' 

/ \ 7-2 

dt ^ dx m, dx m.n.n \ mn I dx ^ V / 



dt ^ dx rrii dx miUio \ n-ig J dx ' 

Hie (dv, I ^ dv^\ ^ J_dl ^dPe I _^_d_ ( _i_&\/nl\ 

rrii \ dt ^ dx J nii dx mine dx 2memi dx \ ^Jn^ dx^ J 



^ = Aneirie - rii), 

In order to obtain the normalized set of equations, following scalings can be used 

c t —Tn 

X — )■ ^-x, t — )■ — , n — )■ nno, u — )■ mCs, — )■ — - — . (2) 

UJpi UJpi 6 



where, ujpi = ^jATre^neo/mi and Cs = \fm^^(?Jrni are characteristic plasma frequency and 
sound-speed, respectively. Also, rigo denotes the electron equilibrium density. Furthermore, 
the electron degeneracy pressure in fully degenerate configuration can be expressed in the 
following form [l2l | 
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|i?(2i?2-3) Vl + i?2 + 31n i? + Vl + i?2 |, (3) 

in which R = ppe/^eC = (^e/'^o)^^'^ (^0 = ^"^^3'^ — 5.9 x lO^^cm"^), where, ppe is the 
electron Fermi relativistic-momentum. It is noted that in the limits of very small and very 
large values of the relativity parameter, R, we obtain (e.g. see Ref. 12|) 

20 U;* me""' 

HI)' hen 




Pe={ '^^'^i '"-^r ■ }. (4) 



Therefore, the relativistic normahzed hydrodynamics set of plasma equations including the 



electron tunneling term, considering the fact that l/ni.{dPe/ dx) = dyl + B? / dx, may be 
written as 

duj _|_ dnjVi _ n. 
dt dx 

at dx 



m&. (dve I y dve\ ^ 91 _ ^^/l+RlT 
rrii \ dt ^ dx J dx dx 



at 'ox ax \ jiinTieC / \ riin i dx ' \ J 



2 V rrii j dx \ ^Jn^ 8x^ 

= n.-n., 



where, we have used R = Ron^J^ and Rq = {neo/noY^^ is a measure of the relativistic 
effects (i.e. relativistic degeneracy parameter) and = h^uj^J {nieC^Y , is the relativistic 

quantum diffraction parameter which is related to the relativistic degeneracy parameter via 
the relation {me/mi)H^ /2 ~ 4.6 x IO^^^-Rq {upi = ^/Aire'^nQRl/mi). It should be noted 
that the parameter Rq is related also to the mass-density (of white dwarf, for instance) 
through the relation p ^ 2mp?2eo or p{gr/av?) = PqRq with pQ{gr/cim?) ^ 1.97 x 10^, 
where, rrtp is the proton mass. The density po is exactly in the range of a mass-density 
of a typical white dwarf. The density of typical white dwarfs can be in the range 10^ < 
p{gr/av?) < 10^, which results in values of 0.37 < -Rq < 8 for the relativity parameter. 
S. Chandrasekhar, in ground of Fermi-Dirac statistics, has proven that for a white dwarf 
with a mass density p, the electron degeneracy pressure turns from oc p^^^ dependence 
for normal degenerac y ( Ro <^ 1) to Pg oc p"^^^ dependence for ultra-relativistic electron 
degeneracy (-Rq ^ 1) ll3|- For the highest value of -Rq = 8 (as for the typical white-dwarf), 
we get the negligible value of [me/mi)!!^ /2 ~ 10^^'^, which validates the applicability of 
current model to white dwarfs or even neutron stars. In the forthcoming analysis higher 
values of the relativistic degeneracy parameter has also been considered, since, for neutron 
stars the electron densities of order Ue — 10^^ {cm~^) is possible. More recently, A. A. Mamun 



and P. K. Shukla [2l| have considered the nonlinear solitary excitations and double-layers in 
dusty degenerate dense plasmas in ultra-relativistic limit, however, as it is noted from Eqs. 
(j4]), their corresponding limit is only a special case, while, the present model includes the 
whole range of relativistic degeneracy parameter, i?o- Now, neglecting the term containing 
the ratio me/mj and assuming that fc^Tj ^ TrieC^, we arrive at the following simplified set 
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of basic normalized equations 

drij _|_ driiVj _ n 
dt dx ^' 

drie _|_ drieVe n 

at dx 

dt ^« dx ~ dx ' 

_ _ 
dx^ — ''^e ■'^i- 

Solving the basic equations (Eqs. ([6])) for the electron and ion number densities in terms 
of electrostatic potential, at the quasi-neutrality condition, and employing the appropriate 
boundary requirement ( lim rig j = 1 and lim = 0), results in the following energy- density 
relations 

1 



/i- 



ne = Ro^ 



24, 



(7) 



Rl + (I)(2^/TTM + 



3/2 



These equations are somewhat different from the standard energy-density relations for 



non-relativistic semiclassical Thomas- Fermi model [20|, |22|, |23|], as expected. Note that, 
unlike for Thomas-Fermi model for plasma, in Eqs. there is no explicit contribution 
for electron (Fermi) temperatures, since, there is a hidden dependence with the electron 
equilibrium quantum number-densities. It is worth mentioning that, there is a widespread 
misconceptions about the independence of the particles equilibrium number- density and 
their corresponding Fermi-temperatures in many recent quantum plasma literatures 24l |. 
which is avoided in current model. The state of a completely degenerate quantum plasma 
is determined through only one of the mentioned quantities since the Fermi-quantities are 
fixed via the particle equilibrium density. In the proceeding section we aim at finding 
the appropriate pseudo-potential which describes the propagation of nonlinear ion-acoustic 
waves in the plasma described by Eq. ([6]). 



III. ARBITRARY-AMPLITUDE NONLINEAR PROPAGATIONS 

To obtain the appropriate Sagdeev pseudo-potentials which describes the possibility of 
ISWs as well as IPWs in a relativistically dense electron-ion plasma, we may reduce to 
stationary frame by moving into the new coordinate ^ = x — Mt, where, M is the Mach 
number being a measure of the wave-speed relative to that of plasma sound, Cg. Substituting 
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Eqs. d?]) into Poisson's equation, multiplying by ^ and integrating once with suitable 
boundary conditions, mentioned in the previous section, we derive the well-known energy 
integral 



where, the desired Sagdeev pseudo-potentials f/(0) read as 

= I !^8M'Rl + Ro {2Rl - 3) y^T+Rl - SM^R^^l-^ 



(8) 



+ 3sinh"^i?o + In i 



-3 In 2 



Note that the values of f/(0) and d^U{(l)) both vanish at 
we may also derive the maximum values of as 

M2 



(9) 



0, as required. From Eqs. ([7]) 



(10) 



0^2 = ±l-^l + i?2, (11) 

The possibility of solitary waves then requires the following conditions to met, simultaneously 



^(0)1^=0 



0. 



c/2f/(0) 



4>=o 



< 0. 



(12) 



4>=o 



It is further required that for at least one either maximum or minimum nonzero 0-value, 
we have ?7(0m) = 0, so that for every value of in the range 0^ > > or 0^ < < 0, 
U (0) is negative. In such a condition we can obtain a potential minimum which describes 
the possibility of a ISW propagation. 

On the other hand, the existence of a nonlinear periodic wave requires that 



f^(0)L=o 



c/t/(0) 



0. 



rf2f/(0) 



</>=0 



>0, 



(13) 



</>=o 



It is also required that for at least one either maximum or minimum nonzero 0-value, we 
have f/(0m) = 0, so that for every value of in the range 0m > > or 0^ < < 0, 
f/(0) is positive. In such a condition we can obtain a potential minimum which describes 
the possibihty of a IPW propagation. 
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Evaluation of the second derivative of the Sagdeev pseudo-potential immediately yields 
a critical Mach-number value of 

M„ = "o^^+^r"" , (14) 

below/above which lASW/IAPW may exist. This value can also be found considering the 
nature of these nonlinear waves whose amplitude tends to zero as the Mach-number M tends 
to its critical value, M^j. which may be predicted by expanding the pseudo-potential U (0) to 
the third order in a Taylor series in 5(f) near = 0. However, since, the other conditions on 

the pseudo-potential requires analytic solutions, which is nontrivial to obtain in our case, in 
the proceeding section a numerical scheme should be used. 



IV. NUMERICAL ANALYSIS AND DISCUSSION 

Figure 1 shows the region of existence for ISWs/IPWs in yellow/blue (right/left region) 
in two different scales in M-Rq plane. It is remarked that the allowed Mach-range is signif- 
icantly affected by the relativity parameter, Rq. It is also observed that the IPWs posses 
relatively lower Mach-values compared to ISWs for all values of the relativistic degener- 
acy parameter, Rq. It is further observed that as the relativity parameter increases the 
Mach-range for both nonlinear excitation (ISWs and IPWs) increase and the Mach-number 
range tend to higher values. Figure 1(a) shows that supersonic solitary excitations occur 
beyond the values Rq > 1.3. This value {Rq > 1.3) corresponds to a mass-density of 
p ~ 4.3 X 10^{gr/cm^) which is within the range of mass-density of typical white dwarfs 
(10^ < p{gr/cm?) < 10^). It should be noted, however, that the parameter M = Vg/cg 
is the relativistic Mach-number, since, the normalization parameter Cs ~ 8 x 10^(cm/s) 
is comparable to the Fermi-speed of electrons in typical metals and is much higher than 
the ion-acoustic speed in the non-relativistic non-degenerate plasmas. It should be noted 
also that for all the values of the relativistic degeneracy parameter, Rq, there is a lower 
Mach-number hmit for the propagation of both types of nonhnear excitations (ISWs and 
IPWs). On the other hand. Fig. 1(b) shows that both nonlinear excitations lie entirely in 
supersonic range for relativistic parameters i?o > 25. This value {Rq > 25) corresponds to 
a mass-density of p ~ 3 x 10^° (gfr/cm^), which is larger than that for even a massive white 
dwarf. 
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The profiles of pseudo-potentials and their variation with respect to the relativity pa- 
rameter and Mach-number are shown in Fig. 2 for both ISWs and IPWs. It is remarked 
that, the increase of the relativistic degeneracy parameter Rq decreases/increases the po- 
tential depth and width for lASWs/lAPWs. Also, it is observed that the increase in the 
Mach-number increases/decreases the potential width and depth for ISWs/lPWs. On the 
other hand. Fig. 3(a) shows the solitary excitation profiles for different relativistic degener- 
acy parameter values for a fixed Mach-number. It is confirmed that the amphtude of ISW 
reduces with increase of Rq (which is directly related to the mass-density) , while it increases 
with increases in the value of the Mach-number for a fixed Rq value (Fig. 1(b)). Therefore, 
it is confirmed that the narrower /wider the solitary wave the taller/smaller it is, which is a 
general nonlinear wave feature in agreement with the relation 0qA = est., where, 0o and A 
denote the amplitude and the width of nonlinear waves, respectively. 



V. CONCLUDING REMARKS 



The Sagdeev pseudo-potential approach was used to investigate the propagation of nonlin- 
ear ion waves in a superdense electron-ion plasma with relativistically degenerate electrons. 
The matching criteria of existence of such solitary excitations were numerically investigated 
in terms of the relativity measure of electrons and the allowed Mach-number range for 
propagation of such waves was evaluated numerically. It was shown that the relativistic 
degeneracy of electron has crucial effects on nonlinear wave dynamics in relativistically de- 
generate plasmas such as that found in white dwarfs, neutron stars and the cores of massive 
planets. 



8 



[1] H. Haug, S. W. Koch, "Quantum theory of the optical and electronic properties of semicon- 
ductors" , World Scientific, 2004, 
[2] S. Son and N.J. Fisch, Phys. Lett. A, 329, 76(2004). 
[3] S. Son and N.J. Fisch, Phys. Rev. Lett., 95, 225002(2005). 
[4] R. G. Greeves and C. M. Surko, Phys. Rev. Lett. 75 3846(1995). 

[5] Shapiro, S. L. and Teukolsky, S. A., Black Holes, White Dwarfs and Neutron Stars, (Wiley, 

New York), (1983). 

[6] M. J. Rees, in: G. B. Gibbons, S. W. Hawking, S. Siklas (Eds.), The Very Early Universe, 

Cambridge Univ. Press, Cambridge, (1983). 
[7] H. R. Miller, P. J. Witta, Active Galetic Nuclei, Springer- Verlag, Berlin, p.202, (1987). 
[8] P. Goldreich, W. H. Juhan, Astrophys. J., 157, 869(1969). 
[9] F. C. Michel, Rev. Mod. Phys. 54, 1(1982). 
[10] M. Bonitz, D. Semkat, A. Filinov, V. Golubnychyi, D. Kremp, D.O. Gericke, M.S. Murillo, 

V. Filinov, V. Fortov, W. Hoyer, J. Phys. A, 36 5921(2003). 
[11] Shukla P K, Eliasson B, "Nonlinear aspects of quantum plasma physics" Phys. Usp. 51 
53(2010). 

[12] S. Chandrasekhar, "An Introduction to the Study of Stellar Structure", Chicago, 111. , (The 

University of Chicago press), (1939) pgs. 360-362. 
[13] S. Chandrasekhare, Science, 226 4674(1984). 

[14] W. Misner, K.S. Throne, J.A. Wheeler, Gravitation, Freeman, San Francisco, p.763, (1973). 
[15] R. Sabry, W. M. Moslem, F. Haas, S. AH, and P. K. Shukla, Phys. Plasmas, 15, 122308(2008). 
[16] A. P. Misra, N. K. Ghosh and C. Bhowmik, Eur. Phys. J. D, 49 (2008)373 
[17] A. P. Misra, P. K. Shukla, C. Bhowmik, Phys. Plasmas, 14, 082309(2007). 
[18] Prasanta Chatterjee, Kaushik Roy, S. V. Muniandy, and C. S. Wong, Phys. Plasmas 16, 
112106(2009). 

[19] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 052302(2010). 

[20] U. M. Abdelsalam, W. M. Muslcm and P. K. Shukla, Phys. Lett. A, 372 4057(2008). 
[21] A. A. Mamun and P. K. Shukla, Phys. Lett. A, 374, 4238(2010). 
[22] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 072101(2010). 



9 



[23] M. Akbari-Moghanjoughi, Phys. Plasmas, "Distinctive features of ion-acoustic solitons in EPI 
super-dense magneto-plasmas with degenerate electrons and positrons" 17 (in-press due Sep. 
2010). 

[24] M. Akbari-Moghanjoughi, "Comment on article 'The effects of Bohm potential on ion-acoustic 
solitary waves interaction in a nonplanar quantum plasma' " [Phys. Plasmas 17, due Sept. 
(2010)] 



10 



FIGURE CAPTIONS 

Figure- 1 

(Color online) The stability regions of arbitrary-amplitude lASWs/IAPWs is shown in 
M-Rq plane for a relativistic degenerate electron-ion plasma. The lower region (yellow) 
correspond to the stabihty of ISWs and the upper (blue) region to IPWs. 

Figure-2 

(Color online) The variation of Sagdeev pseudo-potential with respect to different rela- 
tivity parameter, Rq, and Mach- number, M for both ISWs and IPWs. 

Figure-3 

(Color online) The profiles of solitary waves and the variation of their shape with respect 
to different relativity parameter, Rq, and Mach-number, M. 
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